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ABSTRACT 

The absolute luminosities and homogeneity of early-time infrared (IR) light curves of 
type la supernovae are examined. Eight supernovae are considered. These are selected 
to have accurately known epochs of maximum blue light as well as having reliable 
distance estimates and/or good light curve coverage. Two approaches to extinction 
correction are considered. Owing to the low extinction in the IR, the differences in 
the corrections via the two methods are small. Absolute magnitude light curves in the 
J, H and if-bands are derived. Six of the events, including five established "Branch- 
normal" supernovae show similar coeval magnitudes. Two of these, SNe 1989B and 
1998bu, were observed near maximum infrared light. This occurs about 5 days before 
maximum blue light. Absolute peak magnitudes of about -19.0, -18.7 and -18.8 in 
J, H K respectively were obtained. The two spectroscopically peculiar supernovae 
in the sample, SNe 1986G and 1991T, also show atypical IR behaviour. The light 
curves of the six similar supernovae can be represented fairly consistently with a 
single light curve in each of the three bands. In all three IR bands the dispersion 
in absolute magnitude is about 0.15 mag, and this can be accounted for within the 
uncertainties of the individual light curves. No significant variation of absolute IR 
magnitude with i?-band light curve decline rate, Ami5(B), is seen over the range 
0.87 < Ar7ii5(B) < 1.31. However, the data are insufficient to allow us to decide 
whether or not the decline rate relation is weaker in the IR than in the optical region. 
IR light curves of type la supernovae should eventually provide cosmological distance 
estimates which are of equal or even superior quality to those obtained in optical 
studies. 
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1 INTRODUCTION 
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It is increasingly believed that the optical light curves of 
type la supernovae (SNe la) differ between events in intrin- 
sic peak luminosity and shape. Moreover, the SN la peak 
luminosity appears to correlate with its optical decline at 
early times; slower declining SNe la have greater peak lu- 
minosity (Pskovskii 1977, 1984; Phillips 1993; Hamuy et al. 
1995; Riess, Press & Kirshner 1995; Hamuy et al. 1996a; 
Riess et al. 1998; Saha et al. 1999; Phillips et al. 1999). 
Work by Hamuy et al. (1996a), Tripp (1998), Saha et al. 
(1999) and Phillips et al. (1999) indicate slopes of between 
0.5 and 0.8 for the variation of the absolute B magnitude 
at maximum light, Msmax, with Ami5(_B), the decline in 
B magnitude at 15 days post-tBmax- {tBmax is the epoch 
of maximum light in the _B-band. Supernova epochs in this 
paper will usually be relative to this fiducial time.) It is sug- 
gested by some that the relation also contains a non-linear 
component (Riess et al. 1998, Saha et al. 1999, Phillips et 



It has been established that, even after correction has 
been made for the decline-rate relations, the absolute mag- 
nitude at tBmax varies with colour. Phillips et al. (1999) 
endeavour to derive decline-rate relations free of the effects 
of host galaxy dust reddening. They find that, in a sam- 
ple which includes events which are noteably red at maxi- 
mum viz. (Bmax — Vmax) > 0.2, a plot of MBmax against 
{Bmax — Vmax) (both decline-rate corrected) yields a slope 
of 3.5±0.4 which is consistent with the Galactic reddening 
law. Even if the sample is restricted to those SNe la hav- 
ing {Bmax — Vmax) < 0.2 (the great majority), Phillips 
(private communication) still finds a slope of 3.9 ± 1.3. This 
tends to support their contention that most of the decline- 
rate corrected colour variation is due to dust reddening. 

Tripp (1998) and Saha et al. (1999) examine samples 
restricted to {Bmax — Vmax) < 0.2. These are very simi- 
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lar to the {Bmax — Vmax) < 0.2 subset of the PhiUips et 
al. sample. Unhke Philhps et al., neither Tripp nor Saha et 
al. try to separate out possible competing effects of intrin- 
sic colour variation and dust reddening in these samples. 
For their decline-rate corrected Msmax versus B — V slope, 
Tripp (1998) finds 2.4 ± 0.4 and Saha et al. (1999) find 1.7, 
smaller than the value found by Phillips et al. (although, as 
indicated above, the uncertainties are large). If extinction is 
negligible in the [Bmax — Vmax) < 0.2 SN samples then 
the Tripp and Saha et al. results may indicate that signifi- 
cant intrinsic colour variation in SNe la is present even after 
decline-rate correction. However, if the colour variation has 
a significant component due to dust extinction effects then 
it is less clear what the results of Tripp and Saha et al. sig- 
nify. 

In order to exploit fully the SN la phenomenon as a 
valuable distance indicator, it is desirable to find as many 
ways as possible for dealing with the effects of extinction. In 
this paper I examine light curves in the near-infrared (JHK) 
where it might be hoped that the effects of interstellar dust 
are less. This approach was anticipated in the pioneering 
paper of Elias et al. (1981), where high quality JHK light 
curves were presented for the type la supernovae 1980N, 
1981B and 1981D. Two of these, SNe 1980N and 1981D oc- 
curred in the same galaxy, NGC 1316 (Fornax A). By simply 
shifting the SN 1981B light curves by 0.4 mag, they found 
that the resulting scatter in the individual light curves ap- 
peared to be 0.1 mag or less, which was barely larger than 
the uncertainties in the photometry. This was a lower scat- 
ter than that seen in optical light curves at that time, and 
Elias et al. suggested that this was due to the much reduced 
effects of internal reddening in the IR. They suggested that, 
consequently, IR light curves of supernovae might be useful 
as distance indicators out to a few tens of Mpc. In Elias et 
al. (1985), fiducial IR light curves ("templates") were pro- 
duced, made up of a series of straight line segments. These 
were based on JHK data points for the three supernovae 
from Elias et al. (1981) together with additional IR pho- 
tometry from eight other supernovae. The light curves of 
individual supernovae were adjusted in relative magnitude 
to provide the best fits to the fiducial curves. Relative epochs 
were fixed by the estimated tsmax. In spite of the sometimes 
large uncertainty in tsmax, they found that the uncertainty 
in the fiducial curves was roughly ±0.03 mag between days 
t=+5 and t=+40. They suggested that SNe la have a dis- 
persion in absolute IR magnitude of ±0.2 mag, and possibly 
±0.1 mag, making them potentially valuable for distance 
determination within the Local Supercluster. 

Before the work of Elias et al., Kirshner et al. (1973) 
acquired the first ever IR measurements of a SN la viz. 
SN 1972E, but only a part of their data are of quality 
comparable to that of modern photometry. Between the 
publications of Elias et al. (1981, 1985) and the advent of 
SN 1998bu, very little IR photometry of SNe la was re- 
ported. While good coverage of the SN 1986G IR light curves 
was achieved by Frogel et al. (1987), this was a highly atypi- 
cal, subluminous event. A few JHK points were obtained for 
SNe 1984A (Graham et al. 1988); 1989B (Kidger et al. 1989; 
Wells et al. 1994) and 1991T (Menzies & Carter 1991; Har- 
rison & Stringfellow 1991). SN 1998bu in NGC 3368 (M96) 



was the first normal type la since the 3 events described 
in Elias et al. (1981) for which reasonably good IR pho- 
tometric coverage was achieved (Jha et al. 1999, Meikle & 
Hernandez (1999), Hernandez et al. 2000). In addition to the 
occurrence of SN 1998bu, an important development since 
the work of Elias et al. has been the establishment of more 
accurate distances to some of the host galaxies. Given these 
developments, it is timely to carry out a re-examination of 
the extent to which SNe la are homogeneous in the IR. 



2 THE SAMPLE 

I examine the IR [JHK) luminosities for eight SNe la viz. 
SNe 1972E, 1980N, 1981B, 1986G, 1989B, 1991T, 1998bu 
(Table 1). These were selected to have accurately known 
epochs of maximum blue light as well as having reliable dis- 
tance estimates and/or good light curve coverage. The pho- 
tometric data are listed in the Appendix (Table 3). Only 
measurements for which the quoted error was no greater 
than ±0.2 mag were included. (Owing to its peculiarity and 
very sparse IR coverage an exception to this rule was made 
for the peculiar SN la SN 1991T.) Elias et al. (1985) also 
present IR photometry for the type la SNe 19711, 71L, 83G, 
83R, 83U, & 84A, as do Graham et al. (1988) for SN 1984A, 
but these were excluded on the basis of the selection criteria. 
In view of the rather unusual spectral energy distribution of 
SNe la in the JHK bands, I make no attempt to correct for 
differences in the detectors, filters or standards used between 
different observatories. 

2.1 Distances 

The adopted distance moduli are given in Table 1, col. 4. 
For SNe 1972E, 1981B, 1989B & 1998bu, HST-Cepheid dis- 
tances to the host galaxies are available. For all except 
SN 1972E I use the recalibrated Cepheid distance moduli 
given in Gibson et al. (2000). I show only their random er- 
rors in Table 1. They also estimate a systematic error of 
0.16 mag, which would affect all their distance moduli by 
the same amount. The distance modulus to NGC 5253, the 
host galaxy of SN 1972E, is still subject to quite substan- 
tial uncertainty, with Saha et al. (1999) and Gibson et al. 
(2000) disagreeing by 0.5 mag. In view of this discrepancy, I 
have adopted the mean of their values, viz. 27.86 ± 0.25 for 
SN 1972E. 

NGC 4527, the host galaxy of SN 1991T, lies close to 
NGC 4536 and NGC 4496, both of which have HST Cepheid 
distances (Fisher et al. (1999) and references therein). I 
therefore adopt the mean distance modulus of these two 
galaxies for NGC 4527. 

SNe 1980N and 1981D both occurred in the same 
galaxy, NGC 1316 (Fornax A). Unfortunately there is not 
yet a Cepheid distance to this galaxy, and the value of the 
distance is currently subject to some disagreement. Ciar- 
duUo, Jacoby & Tonry (1993) use the planetary nebula lu- 
minosity function (PNLF) distance indicator to obtain a dis- 
tance modulus of 31.19 ± 0.07 based on a distance modulus 
of 24.43 for M31. A similar value for the Fornax Cluster 
is obtained by Bureau, Mould & Staveley-Smith (1996) who 
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derive 31.10±0.45 based on an 7-band TuUy-Fisher distance 
relative to the Virgo Cluster, and the HST-Cepheid distance 
to Virgo (Freedman et al. 1994). Madore et al. (1999) have 
determined an HST-Cepheid distance to the Fornax Clus- 
ter galaxy, NGC 1365, finding a modulus of 31.35 ± 0.07. 
They suggest that the distance to NGC 1316 is about the 
same. However, Saha et al. (1999) state that there are rea- 
sons to suspect that NGC 1365 is in the foreground of the 
Fornax Cluster. Ferrarese et al. (2000) use their Cepheid- 
based calibration of the /-band surface brightness fluctu- 
ation (SBF) distance indicator to obtain a distance mod- 
ulus to NGC 1316 of 31.71 ± 0.19 and a Fornax Cluster 
distance modulus of 31.59 ± 0.04. They find that where I- 
band SBF and direct Cepheid distance measurements to a 
galaxy are available good agreement is obtained, while the 
PNLF distance indicator produces systematically smaller 
distances. (Ferrarese et al. also suggest that NGC 1365 is 
towards the front of the Fornax Cluster, similar to the con- 
clusions of Saha et al). Of course, one might use SN 1980N 
itself to estimate its distance. Hamuy et al. (1996b) obtain 
31.60 ±0.15. Equations (17) & (18) in Phillips et al. (1999) 
yield 31.70 ± 0.08, while Saha et al. (1999) quote a value of 
31.84 ± 0.21. However, in a study such as this, where we are 
attempting to investigate the homogeneity of SNe la, it is in- 
advisable to use distances determined from the phenomena 
we are examining. In view of the uncertainty in the distance 
to NGC 1316, I have adopted the Ferrarese et al. (2000) For- 
nax Cluster distance modulus of 31.59 for this galaxy. How- 
ever, given the typical random errors given by Ferrarese et 
al. and the uncertainty of the location of NGC 1316 within 
the Fornax Cluster, I assume a somewhat larger error of 
±0.1 mag in the distance modulus of NGC 1316. 

There are no Cepheid distance measurements for 
NGC 5128 (Centaurus A), the host galaxy of SN 1986G. 
Therefore I used the weighted mean distance obtained from 
Cepheid-calibrated tip of the red giant branch and /-band 
SBF methods (Ferrarese et al. 2000). This gives a distance 
modulus of 28.01 ± 0.12 



2.2 Extinction Correction 

As indicated above, for most SNe la there is disagreement 
as to how much of the decline-rate-corrected colour distri- 
bution can be attributed to dust extinction. I have therefore 
investigated two approaches to correcting the extinction. In 
the first approach I have followed the work of Phillips et al. 
(1999) which provides estimates of E{B — V) for seven of the 
supernovae in my sample. They did not consider SN 1981D, 
but Mark Phillips kindly provided his estimate for this su- 
pernova, using the same technique as in Phillips et al. (1999). 
He obtains a total E{B-V) = 0.18±0.09. The IR data were 
then de-reddened assuming the extinction law of Cardelli, 
Clayton & Mathis (1989) with Rv = 3.1. 

In the second approach, no reddening corrections were 
made to the four supernovae for which {Bmax — Vmax) < 
0.2. For the four that exceeded this value (SNe 1981D, 
1986G, 1989B & 1998bu) I made a correction for extinc- 
tion using the Cardelli, Clayton & Mathis law. However, I 
used a lower value of Rv, viz. 2.6, which is the value found 



by Phillips et al. (1999) from a fit to colour-magnitude plots 
for events for which they derived non-zero E{B — V). The 
estimated total extinctions are shown in Table 1, cols. 6-11. 



3 LIGHT CURVES AND COLOUR 
EVOLUTION 

In this section, I examine the individual light curves and 
colours of the eight supernovae. For each supernova I con- 
verted the photometry to intrinsic absolute magnitudes us- 
ing the distance moduli and extinctions given in Table 1. It 
was found that the two approaches to dealing with the ex- 
tinction correction made little qualitative difference to the 
conclusions. Of course, a major reason for this IR study was 
to try to reduce effects of extinction uncertainties. As would 
be expected, the luminosities derived by method (1) were 
slightly higher. Further details are given in Section 4. In the 
presentation of the absolute light curves and colour evolu- 
tion that follows, de-reddening is by means of method (2). 

3.1 Light Curves 

In Figures 1, 3 & 5 I show the intrinsic absolute magnitude 
IR light curves for the J, H and K bands covering the pe- 
riod up to +110 days. In order to display clearly the much 
greater density of data at early times, I show the same light 
curves in Figures 2, 4 & 6, but plotted out to just +60 days. 
Also plotted are the template IR light curves of Elias et al. 
(1985). These are based on JHK light curves of SNe 1972E, 
80N, 81B & 81D plus a few points from other SNe la. I 
adjusted the templates in epoch to give the best match to 
SN 1980N. This required Elias et al.'s fiducial time, "to = 0", 
to be set to -6.25 days (i.e. 6.25 days before tsmax). I also 
truncated the templates at +2 days since this is the earli- 
est epoch for which Elias et al. (1981, 1985) presented data. 
The vertical positions of the templates are set at the aver- 
age absolute magnitude values for six of the events. This is 
discussed in detail in Section 4. It can be seen that, with the 
exception of SN 1986G, the templates provide a useful rep- 
resentation of the IR light curves and that the dispersion in 
coeval magnitudes is not great. I now examine the IR light 
curves by era. 

3.1.1 t < tsmax 

There are four SNe for which data in this era are available 
VIZ. SNe 1986G, 1989B, 1991T & 1998bu. The earliest-ever 
IR measurement of a type la supernova is the single JHK 
observation of SN 1991T by Menzies & Carter (1991) at 
about -11 days. This is discussed later. 

The observations of SN 1998bu show that SNe la peak 
in the IR about 5 days before tsmax. At t = — 5 days we 
have data for SNe 1986G, 1989B and 1998bu. SNe 1989B 
and 1998bu have a similar J-magnitude of about -19.0 at - 
5 days. At H &c K, SN1998bu peaks at about -18.7 in H and 
-18.8 in K, while SN 1989B is brighter by about 0.25 mag. 
Given the uncertainty in the distances of SNe 1989B and 
1998bu, this difference is not significant. SN 1986G is 0.4- 
0.5 mag fainter than SN 1998bu. Random errors in distance 
moduli and extinction corrections are probably insufficient 
to account for the apparent lower luminosity of SN 1986G. 
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I conclude that this supernova was indeed significantly un- 
derluminous at this time. 

3.1.2 tBmax < t <+25 d 

Since the work of Elias et al. (1981) it has been known that 
Type la supernovae light curves exhibit two maxima in the 
JHK bands. A similar double peak is also seen in the /-band 
(Ford et al. 1993). In the J-band, a particularly pronounced 
mininum is seen at t = +15 d (c/. Figs. 1 & 2). In the 
H and K-ha.nAs (Figs. 3-6), a less pronounced minimum 
is seen, occurring somewhat earlier, at about t = +10 d. 
Interestingly, in spite of its peculiar nature, the light curves 
of SN 1986G exhibit absolute magnitudes similar to those 
of the other supernovae around +7 days. Nothing can be 
said about SNe 1972E or 1991T as no IR observations are 
available during this time. 

3.1.3 +25d<t<110d 

For three of the four events for which we have post-+25 d 
light curves, we can see that the decline rates and coeval 
absolute magnitudes are similar. The exception is SN 1986G 
which is fainter in JHK by typically one magnitude during 
much of this period. It is unlikely that such a large factor is 
the result of errors in the estimates of distance or extinction. 

3.2 Colours 

In Figs. 7 & 8 I plot the (J - H) and {H - K) evolution, 
respectively. As above, I consider only the case in which the 
extinction was corrected by method (2). Also in Figs. 7 & 8 
I show the colour evolution templates of Elias et al. (1985), 
again truncated at day +2. 

3.2.1 {J -H) colour 

The t=-ll d point of SN 1991T is at about {J - H) = +0.6. 
However, given that this is the only point at this epoch and 
that 91T had an unusual optical spectrum at early times 
(Fisher et al. 1999), such a red colour may not be typical. 
Between -8 days and +2 days, data exist for only SNe 1986G 
and 1998bu, plus SN 1989B at a single epoch. There is little 
evidence of colour change in this period. SN 1998bu remains 

at {J ~ H) ' 0.2 with SN 1986G somewhat redder at 

~ —0.05. (It should be remembered that these colours have 
already been corrected for extinction). Between +2 d and 
+ 10 d the behaviour of the six supernovae observed is quite 
uniform. At +4 d the (J — //) colour moves sharply to the 
red, reaching {J - H) ^ +1.0 by day +10. Even SN 1986G 
is not strikingly different in its colour evolution during this 
period. Between +10 and +15 days, the general trend is 
continued reddening, reaching (J — H) — +1.3 by day +15. 
The colour then moves to the blue reaching +0.65 by day 30. 
After that the colour reddens again, levelling out at about 
+ 1.6 by +65 d. The post-+2 d evolution is essentially that 
described by Elias et al. (1985). The only clear exception to 
this behaviour is that exhibited by SN 1986G. It begins to 
move back to the blue at only +10 days, and then to the 
red again at only +20 days. SN 1986G seems to mimic the 
(,/ — H) behaviour of the other supernovae, but at a higher 



rate of evolution, and never attaining quite the same amount 
of redness. 

3.2.2 (H-K) colour 

Again, the very early SN 1991T point shows an apparently 
unusually red colour. During the 8 days before tsmax, the 
{H — K) evolution of SN 1998bu shows a gradual redden- 
ing from -0.1 to +0.25. After tBmax, there is a gradual 
move to the blue, reaching {H — K) = —0.15 at +40 days. 
The behaviour of SN 1986G is again different. Up to about 
+50 days it is slightly bluer than typical. From the very few 
post-+80 d observations, we note that SNe 1980N & 1981B 
show a gradual drift to the red again. In contrast, SN 1986G 
shows a very striking move to the blue. The unusually blue 
{H — K) behaviour of SN 1986G was pointed out by Frogel 
et al. (1987). The post-+2 A {H ~ K) behaviour of the other 
supernovae is essentially as described by Elias et al. (1981). 

The non-monotonic behaviour of the IR light curves 
and the development of the IR colours have been consid- 
ered by a number of authors. Hoflich, Khokhlov & Wheeler 
(1995) have found that in some of their models, the second 
maximum in the light curve is reproduced in some regions 
of the IR. They attribute this to a time-dependent opac- 
ity causing the effective emitting area to peak at a delayed 
time in the IR compared with the optical region. However, 
they have difficulty in reproducing the strength of the effect 
as observed in all bands from I to K. From analyses of IR 
spectra, Spyromilio, Pinto & Eastman (1994), and Wheeler 
et al. (1998) conclude that the appearance of the red (J — H) 
colour can be attributed to the relatively large reduction of 
line-blanketing opacity in the 1.2 — 1.5 ^m region. Wheeler 
et al. show that the depth of the J-band deficit can provide 
a valuable temperature diagnostic for the silicon layers. 



4 QUANTITATIVE COMPARISON OF 
ABSOLUTE MAGNITUDES 

The presentation given above suggests that the IR behaviour 
of SNe la is generally quite homogeneous. I now endeavour 
to quantify this behaviour. A problem is that the coverage 
of IR light curves is still much inferior to that achieved in 
the optical region. In particular we cannot, as optical stud- 
ies do, compare directly the absolute magnitudes at maxi- 
mum light. However, even in optical studies maximum light 
is sometimes missed. Nevertheless, the magnitude at max- 
imum light can still be estimated by fitting template light 
curves to the data {e.g. Hamuy et al. 1995). I have adopted a 
similar approach for this study. As pointed out in the previ- 
ous section, the IR light curve templates of Elias et al (1985) 
provide a fair representation of the IR light curve shapes. I 
have therefore used these templates to compare the IR ab- 
solute magnitudes of the supernovae. 

I first compared the _ff-band template with the 
SN 1980N iJ-band photometry and varied the position of 
the template in both axes (absolute magnitude and time) to 
minimise xi- For the 14 points considered, a xt minimum 
of 1.7 was achieved by setting to = at —6.25 ± 1.0 days, 
where to is the epoch as defined by Elias et al. (1995). This 
is consistent with the value of -6 days adopted by Elias et al. 
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I fixed all three templates such that to = is at -6.25 days. 
I then adjusted the vertical (absolute magnitude) position 
of the templates to provide a minimum x^-fit to the data for 
all supernovae and all 3 bands. In spite of the fact that the 
templates were originally largely based on SNe 1972E, SON, 
81B & 81D, some of the best fits to these particular SNe 
had values of xt large as 10. However, it was found that 
modest increases in the quoted errors on the data points (up 
to ±0.1 mag) would bring xt close to unity. For the 10-12 
points of SN 1998bu, the fits were less satisfactory, yielding a 
xi exceeding 10. Satisfactory fits were only obtained by arti- 
ficially increasing the errors to around ±0.15 mag. The need 
to arbitrarily increase the individual errors to achieve satis- 
factory fits probably indicates that the IR light curves are 
not completely homogeneous in shape. For SN 1989B good 
fits were easily obtained for the 2 points which overlapped 
the timespan of the templates. For SN 1986G, in view of its 
clearly different light curve shape in all three bands, no tem- 
plate fits were performed. Instead, the weighted mean values 
for the 4 points in the ±12 to ±15 day period were used to 
estimate its Mi4(J), M\4,{H) and M\4[K) (defined below). 
For SN 1991T only a single epoch [t = ±55 d) coincided with 
the template timespan and so the templates were simply ad- 
justed to coincide with the single point at each waveband. 
However, given the peculiar nature of SN 1991T, it is possi- 
ble that the templates are not representative of its true IR 
light curves. 

In Table 2 I show the results of the template fitting 
procedure. The absolute magnitudes (cols. 2, 4 & 6) are 
expressed in terms of Mi4(J), M\i{H) and M-n{K), the 
absolute magnitudes a.t t = ±13.75 days. This epoch cor- 
responds to the epoch of the H — H20 ~ 0.0 point on the 
original H-band template of Elias et al. (1985). It lies close 
to the epoch of first minimum in the J-band. The tabulated 
values are those obtained after increasing the errors on the 
data points until was about unity. (When the fitting pro- 
cedure was performed with the original quoted errors, the 
absolute magnitudes returned never differed from the values 
shown by more than a few percent.) The numbers in paren- 
theses give the error in the smallest two decimal places of 
the Mi4 values returned by the x^-fitting procedure. The 
only exceptions are SNe 1986G and 1991T where the errors 
in parentheses are based simply on the published photom- 
etry errors (see above). However, given the peculiar nature 
of SN 1991T and the fact that the epoch of observation is 
over a month beyond the fiducial t = ±13.75 day epoch, 
it is possible that the quoted errors on the M14 values are 
significantly underestimated for this supernova. Columns 3, 
5 & 7 give the total random error in the absolute magni- 
tude, comprising the x^-fitting or photometry error (shown 
in parentheses in cols. 2, 4 & 6), random error in the distance 
(c/. Table 1), and the extinction correction error. The error 
in the extinction correction was taken to be half of the dif- 
ference in extinction correction using methods (1) and (2). 

Following the work of Phillips et al. (1999) and others, 
in Figure 9 I compare the intrinsic absolute IR magnitudes 
at t = ±13.75 d, Mi4(J), Mii{H), Mi^K), with AmisiB), 
the decline in the B-band 15 days after maximum blue light 
(Table 2, col. 8). For six of the eight SNe (the exceptions are 
SNe 1986G & 1991T) the AI14 values are reasonably consis- 



tent with zero difference between the events. Single parame- 
ter fits to the Mi4 values for the six events yields xt values of 
1.6, 1.8 and 1.4 in J, H and K respectively. The actual val- 
ues {i.e. the weighted means) are Mi4(J) — —16.86 ± 0.06, 
Mi4.{H) = -18.22 ± 0.05 and MuiK) = -18.23 ± 0.05. 
These values are indicated by the horizontal dotted lines 
in Fig. 9. Henceforth, these six supernovae will be referred 
to as the "IR-normals" . The dispersion in all 3 bands is 
about 0.15 mag for the IR-normals. As a check, I fitted the 
templates to the light curves of all six IR-normals simulta- 
neously (about 60 points per band). The same M14 values 
were obtained. SN 1986G has similar Mi4{H) and Mi4,{K) 
to those of the IR-normals, but is is about 0.5 mag brighter 
in J than the IR-normal mean. However, Figs. 1 to 6 show 
that SN 1986G was clearly underluminous by a substan- 
tial amount before ~0 days and after ~±25 days. At later 
epochs (post-±25 days), SN 1986G is fainter by as much as 
one magnitude in the three bands. SN 1991T appears to be 
overluminous but at a low significance - ~ 2(t in each band. 
It should also be recalled that SN 1991T has only a single 
set of JHK observations (at t — ±55 d) which lies within 
the template temporal span. 

Similar results are obtained when the entire analysis is 
repeated using method (1) for the extinction correction. The 
mean M14 values for the six IR normals are, respectively, 
0.09, 0.04 and 0.02 mags brighter in JHK, with slightly 
smaller xt values viz. 1.3, 1.5 and 1.2. 



5 DISCUSSION AND CONCLUSION 

For all three IR wavebands the template analysis indi- 
cates that six of the eight supernovae considered have sim- 
ilar coeval magnitudes. Indeed, to within the uncertainties, 
they are indistinguishable. Of the six, SNe 1972E, 1981B, 
1989B, 1998bu and, probably, SN 1980N are all spectroscop- 
ically "Branch-normal" (Branch, Fisher & Nugent 1993). 
SN 1981D cannot be classified in this way due to insuffi- 
cient spectroscopic coverage. In contrast, the two exceptions, 
SNe 1986G and 1991T, have long been recognised as spec- 
troscopically peculiar. Nevertheless, in spite of the peculiar 
behaviour of SN 1986G, there is evidence that at around 
±7 days it had a similar IR luminosity to that of the IR- 
normals. 

The IR-normals span 0.87 < Ami5(B) < 1.31. Over 
this range the dispersion in absolute magnitude is about 
0.15 in all three bands, and this can be accounted for al- 
most entirely by the uncertainties. In other words, at this 
level of uncertainty, the IR-normals show no systematic vari- 
ation in absolute J, H or K magnitude with Ami5(B). For 
0.87 < Ami5(-B) < 1.31 both Phillips et al. (1999) and Saha 
et al. (1999) give peak absolute magnitude ranges of about 
0.35 in B, falling to about 0.2 in /. It is possible that the 
absolute magnitude range continues to decline as we move 
further into the IR, becoming too small to detect in JHK 
with our limited data. As a check, I examined the B and V 
band absolute peak magnitudes specifically for the six IR- 
normals. For SNe 1972E, 1981B, 1989B & 1998bu I used the 
absolute magnitudes given in Gibson et al.. For SNe 1980N 

6 1981D I used the 7-band SBF Fornax Cluster distance 
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(Ferrarese et al. 2000) together with the peak apparent mag- 
nitudes given in Hamuy et al. (1991). Extinction corrections 
were made following both methods (1) and (2) of section 
2.1. In both the B and V bands I find no significant trend of 
absolute B or V magnitude with Ami5(_B). The dispersion 
is similar to that seen in the IR. It is therefore likely that the 



vide more reliable estimates of cosmological distances. At 
the lowest redshifts dominated by the Hubble expansion, 
{z = 0.01 — 0.1), the peak IR magnitudes of SNe la would 
fie in the range +14 to +20. These magnitudes are well 
wit hin the range of a 4 m class telescope such as U KIRT 



size of :he sample considered here is simply too small and 



the uncertainties in the absolute IR magnitudes too large to 
reveal a Ami5(_B) relation of even comparable strength to 
that found in the optical region. 

The IR magnitudes of SN 1991T at -11 days are J = 
+ 11.73, H = +11.16 and K = +10.93 (Menzies & Carter 
1991). These translate to the highest-ever IR- luminosities for 
a type la event, especially in the H &i K region (c/. Figs. 1- 
6). Although we lack coeval points from other supernovae, it 
appears that SN 1991T was exceptionally IR-luminous for a 
SN la at this epoch. Menzies & Koen (1991) report simulta- 
neous UBV photometry yielding V = +12.4, B-V = +0.11 
and U ~ B = -0.71. Lira et al. (1998) report similar BVR 
magnitudes obtained just 0.3 d after the Menzies & Carter 
IR measurement, giving V — +12.45, B — V = +0.09 and 
V — R — +0.06. Thus, the IR magnitudes indicate a strong 
excess with respect to the optical region e.g. V — K +1.5. 



(c/. |http://www .jach.hawau.edu/JACpublic/UKIRT/ 
instrumcnts/ufti/sensitivities.htm^ ). Indeed, the nearer su- 



pernovae could be successfully monitored by only a 
2 m class telescope such as the Liverpool Telescope (c/. 
http:/ /telescope. livjm.ac.uk/inst /index. html ) . At higher 
redshifts, say 0.5, rest frame J-band and H-band emission 
are shifted into roughly the H- and K-hands respectively. 
(Rest frame K-hand emission would be lost in the thermal 
IR region.) In a flat, A = universe an intrinsic peak ab- 
solute J-magnitude of -19.0 yields an apparent magnitude 
of about +22.3. This estimate includes only an approximate 
K-correction (Poggianti 1997). As more SNe la IR spectra 
become available, it should be possible to make more accu- 
rate estimates of the K-corrections. A J ~ +23 SN la at 
z — 0.44 has been successfully detected by J. Spyromilio 
in a 1-hour integration with the 3.5 m ESO NTT (Riess 
et al., in preparation). Clearly, even deeper coverage would 
be possible with an 8 m class telescope such as Gemini ( cf. 



http : / / www, ast . cam .ac.uk /sciops/instruments /niri /NIRl 
Index. litm| ) 



The ne rt IR observation of SN 1991T was about 10 days 
later. This took the form of a JHK spectrum (Meikle et al. 
1996), and it shows no evidence of excessive IR luminosity. 
The -11 day photometry of SN 1991T was the earliest ever 
observation of a SN la in the IR. The next earliest are the 
IR spectra of SN 1994D (Meikle et al 1996) and photometry 
of SN 1998bu, both at about -8.5 d. They show no signs of 
unusually high luminosity. It may be that the exceptional 
early-time IR luminosity and redness of SN 1991T are fur- 
ther manifestations of its intrinsic peculiarity. However, for 
the present, this result is something of a puzzle. 

The analysis presented here indicates that most SNe la 
are standard IR candles at a level of about 0.15 mag. As 
already mentioned, Elias et al. (1995) suggested that SNe la 
might have a dispersion in IR absolute magnitude as low 
as ±0.1 mag. The present work is consistent with this, al- 
though it also shows there can be exceptional events which 
are much brighter or fainter than the norm. Moreover, al- 
though I find no indication of differences in coeval magni- 
tudes for the six IR-riormals, there is tentative evidence that 
the IR light curve shapes are not exactly identical. 

It should be kept in mind that this study makes use of 
only eight supernovae, and that their light curves are gener- 
ally much less well-sampled than in the optical region. The 
total number of photometry points for all six IR-normals is 
only about 60 in each band within the span of the Elias et al 
templates. It is clearly desirable to pursue the use of SNe la 
IR light curves as distance indicators. As a first step, we need 
to build up a good set of IR light curves for nearby SNe la 
in order to establish the extent to which they are homoge- 
neous. In addition, this will extend the wavelength range 
over which colour information is available and so should 
help to distinguish unambiguously between intrinsic colours 
and dust reddening. Compared with the optical region, the 
lower sensitivity to extinction uncertainties in the IR should 
mean that observations in this wavelength region will pro- 
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Table 1. Supernova Sample 



Supernova t^^ 



(d) 



Host Gal. 



E{B - V) 



Source ' 













(1) 


(2) 


(1) 


(2) 


(1) 


(2) 




1972E 


41449.0(1.0)0 


NGC 5253 


27.86(25) 


0.07(04) 


0.06 


0.0 


0.04 


0.0 


0.02 


0.0 


1 


1980N 


44585.8(0.5) 


NGC 1316 


31.59(10) 


0.07(05) 


0.06 


0.0 


0.04 


0.0 


0.02 


0.0 


2 


1981B 


44672.0(0.2) 


NGC 4536 


30.95(07) 


0.13(03) 


0.11 


0.0 


0.07 


0.0 


0.04 


0.0 


1,2 


1981D 


44679.9(0.5) 


NGC 1316 


31.59(10) 


0.18(06) 


0.16 


0.13 


0.10 


0.08 


0.06 


0.05 


2 


1986G 


46561.5(1.0) 


NGC 5128 


28.01(12) 


0.61(05) 


0.54 


0.45 


0.34 


0.28 


0.22 


0.18 


3 


1989B 


47565.3(1.0) 


NGC 3627 


30.06(17) 


0.37(04) 


0.32 


0.27 


0.20 


0.17 


0.13 


0.11 


4,5 


1991T 


48375.7(0.5) 


NGC 4536 


31.07(13) 


0.16(05) 


0.14 


0.0 


0.09 


0.0 


0.06 


0.0 


6,7 


1998bu 


50953.3(0.5) 


NGC 3368 


30.20(10) 


0.35(03) 


0.31 


0.26 


0.19 


0.16 


0.12 


0.10 


8,9,10 



" Epoch sources - SN 1972E: Leibundgut et al. (1991), SNe 1980N & 1981D: Hamuy et al. (1991), SN 1981B: Shaefer (1995), 
SN 1986G: Phillips et al. (1987), SN 1989B: Wells et al. (1994), SN 1991T: Lira et al. (1998), SN 1998bu: Hernandez et al. (2000) 
Julian Date - 2400000. 

Distance Modulus: the error in parentheses is the random error only - see text. 

E(B - V) for all SNe except SN 1981D is from Phillips et al. (1999) and includes Galactic reddening. E{B - V) for SN 1981D is from 
M. Phillips (private communication). 

The IR extinctions were determined using two methods: (1) & (2) - see text. 
^ IR photometry obtained from: 

1. EUas et al. (1985), 2. Ehas et al. (1981), 3. Progel et al. (1987), 4. Kidger et al. (1989), 5. Wells et al. (1994), 

6. Menzies & Carter (1991), 7. Harrison & Stringfellow (1991), 8. Mayya, Puerari & Kuhn (1998), 9. Jha et al. (1999), 

10. Meikle & Hernandez (1999); Hernandez et al. (2000). 
3 Numbers in parentheses give random error in the smallest two decimal places. 



Table 2. Absolute magnitudes at t = +13.75 days. 



SN 


Mi4(J)g 


AMi4(J)[ 


Mi4{H) 


AMi4{H) 


Mi4{K) 


AMi4{K) 


Ami5(S)| 


1972E 


-17.26(16) 


0.30 


-18.52(04) 


0.25 


-18.64(06) 


0.25 


0.87(10) 


1980N 


-17.02(03) 


0.11 


-18.29(01) 


0.10 


-18.33(01) 


0.10 


1.28(04) 


1981B 


-16.70(04) 


0.10 


-18.03(03) 


0.09 


-18.09(02) 


0.08 


1.10(07) 


1981D 


-16.94(05) 


0.11 


-18.32(03) 


0.10 


-18.24(02) 


0.10 


1.27(09) 


1986G 


-17.39(04) 


0.14 


-18.27(02) 


0.13 


-18.18(02) 


0.13 


1.73(07) 


1989B 


-16.71(04) 


0.18 


-18.03(06) 


0.18 


-18.24(13) 


0.21 


1.31(07) 


1991t| 


-17.47(30) 


0.33 


-18.63(14) 


0.20 


-18.78(25) 


0.28 


0.94(05) 


1998bu 


-16.81(07) 


0.13 


-18.29(03) 


0.11 


-18.27(04) 


0.11 


1.01(05) 



° Numbers in parentheses give the error in the smallest two decimal places of Mi4(J), M\4{H) and M\4{K) returned by the x^-fitting 
procedure. 

' Total random error in the absolute magnitude. It includes x^-fitting error, random error in the distance estimates (c/. Table 1), and 
extinction correction error (see text). 

Values of Ami5(B) taken from Phillips et al. (1999), apart from SN 1981D. For this supernova the value was provided by M. Phillips 
(private communication). Numbers in parentheses give the error in the smallest two decimal places. 

Given the peculiar nature of SN 1991T and the fact that the epoch of observation is over a month beyond the fiducial t = +13.75 day 
epoch, it is possible that the quoted errors on Afi4(J), Mi4{H) and AIi4{K) are significantly underestimated for this supernova. 
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Table 3. Appendix: Infrared photometry of type la supernovae 



Supernova tBmaxid) 



H 



K 



Source 



SN 1972E 



41459.7 
41461.8 
41463.6 
41472.6 
41482.8 
41487.7 
41488.7 
41493.7 
41496.7 
41502.7 
41516.7 
41517.7 
41518.7 



10.37(12) 
10.76(10) 

11.16(13) 
12.58(16) 



9.32(8) 
9.07(8) 

9.76(8) 
9.83(10) 
10.19(10) 
10.39(10) 
10.61(17) 
11.16(16) 
11.17(11) 



9.44(20) 

9.29(16) 

9.35(8) 

9.18(8) 

9.00(16) 

9.95(8) 

9.96(9) 

10.35(10) 

10.24(11) 

11.40(17) 

11.29(8) 

11.03(13) 



Elias et al. (1985) 



SN 1980N 



44590.7 
44591.5 
44592.7 
44596.6 
44618.6 
44621.6 
44624.6 
44631.7 
44653.5 
44657.6 
44667.5 
44676.5 
44677.5 
44685.5 



13.32(5) 
13.50(7) 
13.61(4) 
14.40(4) 
14.15(3) 
14.47(3) 
14.76(4) 
15.42(5) 
16.74(8) 



17.65(16) 
17.69(20) 
17.54(19) 



13.43 
13.46 
13.38 
13.44 
13.47 
13.66 
13.84 
14.18 
15.11 
15.15 
15.74 
16.00 
16.24 
16.50 



1) 



13.23(8) 
13.32(9) 
13.32(5) 
13.38(4) 
13.59(2) 
13.82(2) 
14.01(3) 
14.34(4) 
15.20(11) 



16.21(11) 



Elias et al. (1981, 1985) 



SN 1981B 



44674.67 

44679.0 

44680.73 

44685.7 

44687.7 

44689.0 

44689.8 

44690.8 

44691.7 

44692.8 

44695.8 

44708.7 

44737.7 

44738.5 

44748.5 

44768.5 

44782.5 



12.70(5) 

13.19(7) 

13.75(5) 

14.18(3) 

14.20(3) 

14.29(2) 

14.11(3) 

14.17(2) 

13.85(3) 

13.99(3) 

13.66(3) 

14.01(11) 

16.09(15) 

16.19(6) 

16.65(5) 

17.52(11) 

17.87(13) 



13.00 
12.95 
13.27 
12.98 
12.90 
12.93 
12.84 
12.82 
12.69 
12.73 
12.59 
13.19 
14.53 
14.53 
15.02 
15.80 
16.22 



12.60(5) 
12.81(5) 
12.99(5) 
12.89(2) 
12.90(2) 
12.87(4) 
12.67(2) 
12.81(3) 
12.77(2) 
12.81(4) 
12.68(5) 

14.67(7) 
14.67(3) 
15.03(5) 
15.71(5) 
15.99(13) 



Elias et al. (1981, 1985) 



SN 1981D 



44683.5 
44685.5 
44686.5 
44687.5 
44688.5 
44693.5 
44695.5 



13.38(3) 
13.74(2) 
14.05(2) 
14.27(3) 
14.44(3) 
14.58(4) 
14.49(6) 



13.41 
13.45 
13.56 
13.56 
13.59 
13.28 
13.22 



13.36(4) 
13.43(3) 
13.55(3) 
13.57(3) 
13.58(4) 
13.28(7) 
13.25(6) 



Elias et al. (1981, 1985) 



" JD - 2400000 

^ Numbers in parentheses give the error in the smallest one or two decimal places. 
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Table 3 (continued). Appendix. Infrared photometry of type la supernovae 



Supernova igmax (d) 



H 



K 



Source 





46556.4 


10.10(3) 


9.96(3) 


9.90(3) 


rrogel et al. (1987) 




46557.4 


10.06(3) 


9.97(3) 


9.87(3) 






46562.2 


10.12(3) 


9.99(3) 


9.87(3) 






46563.4 


10.28(3) 


10.10(3) 


10.03(3) 






46565.8 




10.10(3) 


10.02(3) 






46567.3 


10.73(3) 


9.99(3) 


10.03(3) 






46568.3 


10.86(3) 


10.02(3) 


10.05(3) 






46568.4 


10.81(3) 


10.03(3) 


10.08(3) 






46569.3 


10.88(3) 


9.96(3) 


10.01(3) 
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Kidger et al. (1989) 
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Wells et al. (1994) 
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Menzies & Carter (1991) 
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Harrison &; Stringfellow (1991) 
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Jha et al. (1999) 
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Hernandez et al. (2000) 
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11.49(5) 
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11.63(3) 
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11.79(3) 


11.77(3) 






50976.8 


13.08(1) 
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" JD - 2400000 

^ Numbers in parentheses give the error in the smallest one or two decimal places. 
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Figure 1. J-band light curves of type la supernovae. The absolute magnitudes were derived as explained in the text, using method (2) to 
correct for the extinction. The error bars give the photometry errors only. The data set for each supernova is also subject to uncertainty 
in distance and extinction correction. The continuous line is the template light curve of Elias et al. (1985) with their to = set at 
—6.25 days, and with Mi4(J) = —16.86 (see text). 
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Figure 2. Detail of Figure 1. 
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Figure 4. Detail of Figure 3. 
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Figure 5. X-band light curves of type la supernovae. The template light curve has Mi4,{K) = —18.23. Other details as in Fig. 1 caption. 
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Figure 6. Detail of Figure 5. 
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Figure 7. Evolution of the (J — H) colour of type la supernovae. The colours have been de-reddened by method (2) as described in the 
text. The error bars represent random error only. The data set for each supernova is also subject to uncertainty in extinction correction. 
The continuous line is the template light curve of Elias et al. (1985) with their to = set at —6.25 days. 
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Figure 8. Evolution of tlie {H — K) colour of typo la supornovae. Other details as in Fig. 7 caption. 
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Figure 9. Absolute JHK magnitudes of type la supernovae at +13.75 days (M14) plotted against the blue light curve decline rate 
parameter Am\T^[B). The weighted mean values of Mi4(J), M\4{H) and M\4,{K), for the six "IR-normal" supernovae are indicated by 
the horizontal dotted lines. 
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